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Abstract 

As regular conformal blocks describe the N=2 superconformal gauge the- 
ories in four dimensions, irregular conformal conformal blocks are expected to 
reproduce the instanton partition functions of the Argyres-Douglas theories. In 
this paper, we construct matrix models which reproduce the irregular conformal 
conformal blocks of the Liouville theory on sphere, by taking a colliding limits 
of the Penner-type matrix models. The resulting matrix models have not only 
logarithmic terms but also rational terms in the potential. We also discuss their 
relation to the Argyres-Douglas type theories. 
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1 Introduction 

The Liouville conformal blocks on Riemann surfaces are conjectured to be equal to the 
Nekrasov's instanton partition functions of four-dimensional M = 2, SU(2) supercon- 
formal quiver gauge theories [I], which is now referred to as the AGT relation. The 
complex structure moduli of the (punctured) Riemann surface are identified with the 
marginal gauge couplings of the four-dimensional theory, and the external momenta of 
the vertex operators are now encoded in the mass parameters of hypermultiplets. The 
conformal dimensions of the intermediate states are then interpreted as the Coulomb 
branch parameters of the gauge theory. 

After the above conjecture, the notion of "irregular conformal blocks" was given in 
[2] in order to generalize the AGT relation to asymptotically free gauge theories. The 
asymptotically free theories are obtained by taking a scaling limit of mass parameters 
while tuning the gauge couplings so that the physical quantities of the low energy 
theory remain finite. Such a limit corresponds to a colliding limit of two regular vertex 
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operators in the Liouville conformal block, which gives rise to an irregular state, or an 
irregular vertex operator in CFT. The irregular states constructed in [2] are eigenstates 
of L\ or (Li, L2) with non- vanishing eigenvalues. Their generalization to eigenstates of 
(L n , ■ ■ ■ , I/2n) with non-vanishing eigenvalues was recently given in [3] by considering 
a colliding limit of many regular vertex operators (See also [I]). Since such a general 
colliding limit gives Argyres-Douglas theories in the gauge theory side, it was pointed 
out that the conformal blocks with general irregular states inserted should reproduce 
the instanton partition functions of the Argyres-Douglas theories jU [3] . 

On the other hand, the Liouville regular conformal blocks are known to be written 
as the /3-ensemble of matrix models with logarithmic potentials [3] , where the Liouville 
charge is encoded in the /3- deformation parameter of the matrix model. The external 
momenta of the Liouville vertex operators are now interpreted as parameters in the 
matrix model potential and the matrix size N, while the conformal dimensions of 
the intermediate states are encoded in filling fractions, or eigenvalue distributions, of 
the matrix model. This observation is based on some earlier works on the conformal 
symmetry hidden in matrix models j6j El El [9] . From the AGT viewpoint, this type 
of matrix models are expected to reproduce the Nekrasov partition functions of M = 
2, SU(2) superconformal linear quiver gauge theories, and there have been various 

studies in this direction HOI HU H21 H31 HH H51 HH HSl HH [201 Ell 1221 ESI [2H ESl- 

In this paper, we study the matrix model side of the colliding limit of the Liouville 
vertex operators, and derive a series of matrix models which reproduce the general 
irregular conformal blocks of the Liouville theory on sphere. The colliding limit leads 
to not only logarithmic terms but also rational terms in the matrix model potential. 
In fact, we argue that if the matrix model potential is written as a sum of logarithmic 
and/or rational terms then its partition function reproduce a Liouville conformal block 
with insertions of regular and/or irregular vertex operators. Some of the simplest 
cases were already studied in [UJ and expected to reproduce the instanton partition 
functions of 577(2) gauge theories with two and three flavors. On the other hand, in 
this paper, we conjecture that other general matrix models with logarithmic and/or 
rational potentials reproduce the instanton partition functions of Argyres-Douglas 
theories as well as some asymptotically free theories involving the Argyres-Douglas 
theories as their building blocks. 

The rest of this paper is organized as follows. In section [2j we review the colliding 
limit of regular vertex operators and the property of the irregular states. In section 
[3[ we illustrate how the same colliding limit in the matrix model side gives rise to a 
potential with logarithmic and rational terms. We also specify matrix model potentials 
for general irregular conformal blocks. Among the series of matrix models we have 
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constructed in section [3j we study the D 2n type matrix models in detail in section 
[4j We demonstrate that the eigenvalues of (L n , ■ ■ ■ , L 2n ) for the irregular states are 
interpreted as parameters of the matrix model potential while the actions of L k for 
< k < n are now encoded in the filling fractions of the matrix model. In section [4], 
we also illustrate how our matrix models reproduce the small parameter behavior of 
an irregular state of rank 2. In section |5j we discuss the relation between our matrix 
models and Argyres-Douglas type theories, by using the Hitchin system with irregular 
singularities. We show that our matrix models correctly reproduce the Seiberg-Witten 
curve of the corresponding gauge theories. 

2 CFT and irregular singularity 

In [2] , the notion of an "irregular vector" in the representation of the Virasoro algebra 
was given in the context of the AGT relation pQ. It was defined as the following simul- 
taneous eigenstate in the Verma module for a highest weight state of conformal 
dimension A a = a(Q — a): 

L^JW) = Ail/W), Lal/W) = A 2 |J«), L„> 3 |/ (1) ) = 0. (2.1) 

The explicit expression for \I^) was obtained in [26]. The inner product of this state 
is known to be equal to the instanton partition functions of some asymptotically free 
Af = 2, SU(2) gauge theories [2j [27] . The generalization of this irregular vector to 
more simultaneous eigenvalues was given in [3] as follows (See also [1]). First, note 
that the subalgebra 

[L k , L m ] = (k - m)L k+m fork,m>0 (2.2) 

of the Virasoro algebra implies that a simultaneous eigenstate of L k and L m is also an 
eigenstate of L k + m with a vanishing eigenvalue. Therefore, when we consider a vector 
satisfying 

L k \I {n) ) = A k \I (n) ) for n < k <2n (2.3) 

with non-vanishing eigenvalues A*., we find L k \I^) = for k > 2n. Note that the 
action of L k for k < n is not diagonalized by |/( n )). This vector |JW) is called an 
irregular vector of rank n, which depends on the collection of non- vanishing eigenvalues 
A = {A n ,--- ,A 2n }. 

It is noted that the irregular vector may arise when one considers so-called collision 
limit of many primary fields, a certain limiting process which put many fields at the 
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same point. Let us consider a state obtained from (n + 1) primary fields, \R n ) 
E[fc=o ^(Vk'i Afc)|0). This state satisfies the operator product expansion 



T(z)| J R n ) = V (-. r — + — + regular terms) \R n ) (2.4) 

r=0 



\z-y r ) 2 z-y r dy 



where A r = a r (Q — a r ) is the conformal dimension of the primary field ^(y; A r ). 
Taking the limit y r — > 0, a r — > oo while keeping 

n n k 

c = ^a r , c k = Y^ Yl a rY[(-y Sl ) (2.5) 

r=0 r=0 0<s 1 <--<s k <n(s i ^=r) i=l 

finite, one may have 

2n . n-1 „ j 

T{z)\R n )^ (E^ + E ^ + + regUkr termS ) 1 JW) ' (2 ' 6) 

fc=ra fc=0 

where T(z) = E fc W^ +2 and = lim^o.a^oc (IIo<r<.<»(l/r - y s ) 2a ' as \R n )) ■ 

Note here that our Co is denoted by a in [3], and we sometimes write a as Co. The 
constants 

fe 

A fe = (k + l)Qc k - Y c t°k-i (2-7) 

£=0 

are eigenvalues of L k , where we set Q = unless < I < n. On the other hand, L k 



is of a certain operator satisfying the Virasoro algebraic relation (2.2), whose explicit 
expression is given by 

£ k = A k + Y j (£-k)c e -?—. (2.8) 

Thus, the above colliding limit leads to an irregular vector of rank n satisfying 

for 2n + 1 < k, 

L k \I^ n) ) = { A fc |/( n )) for n < k < 2n, (2.9) 
£ fc |/ (n) ) for < k < n - 1. 

Note that (n + 1) eigenvalues A n , • • • , A 2 „ are now encoded in a, ci, ■ • • , c n . The 
coefficients ci, • - • , c n are identified with the eigenvalues of coherent states in the free 
field construction of the irregular states [3]. For n > 2, the quantities a, Ci, ■ ■ ■ , c n 
are not enough to determine the irregular vector |/( ra )) uniquely. In order to fix such 
an ambiguity, it was proposed in [3] that an irregular vector of rank n > 2 can be 
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recursively constructed from irregular vectors with lower ranks. For example, |/^) i 
proposed to be expanded as 



is 



oo 



|/ {2) (ci,c 2 ;a)> 




5>*|/«( Cl ,e*')), 



(2.10) 



where 



Vi = 2(a — a')(Q — a'), v 2 = (« — a') 




Q+2 a ' + 2 a 



) 



(2.11) 



Here [1^ (ci; a)) are so-called "generalized descendants" which are linear combinations 
of vectors obtained by acting L_f c and Ci-derivatives on \I^\ci;a)). In particular, 
\Iq (ci;a)) = |/^(ci;a)) which can be uniquely determined by c\ and a. It follows 



a' in addition to a,c\,C2- The origin of a' can be understood when we note that the 
regular vector |i?2) originally depends on the conformal dimension of the intermediate 
state. 

By generalizing the above argument, we can consider conformal blocks with many 
irregular vertex operators inserted, which are called the "irregular conformal blocks." 

3 Penner model and irregular singularity 

In this section, we derive matrix models which reproduce the irregular conformal blocks 
of the Liouville theory. We start with the matrix model for regular conformal blocks 
[3] and take the same colliding limit in the matrix model side. 

3.1 Penner type matrix models for regular conformal blocks 

Let us first consider the (n + 2)-point correlation function of the Liouville field theory 



which is given by the product of holomorphic and anti-holomorphic regular confor- 
mal blocks. One may evaluate the correlation explicitly in perturbative expansion of 
the Liouville theory using the free correlation ( e 2a iH z ) e 2a 2<t>( w )} = \ z _ y^^aic^ _ rp^g 
correlation is not vanishing if the Liouville charges satisfy the neutrality condition 




from (2.10) that the rank 2 irregular vector depends on an additional parameter 




(3.1) 



n 




(3.2) 



6 



where N is the number of the screening charges so that the correlation is given by 




dXjdXj e 



2b<f>(\ 1 ) 



-,2a a 



(oo) e 2Q *(0) 



n 

k=l 



,2ctu<t>{zk) 



(3.3) 



Coulomb gas 



The holomorphic part of the correlation, or conformal block, is then written as 



■Z/3_Penner J [ [Zk 
0<k<£<n 



-2a k a e 



(3.4) 



up to divergent prefactor, where Zp-p erinei is the partition function of the following 
/3-deformed matrix model [3] 



J /3— Penner 



N 



Y[dX! 



i=i 



A(A) 2/3 exp 



£V(A, 



(3.5) 



Here A(A) = fli 



<KJ<N 



(A/ — Aj) , a//3 = — ib, and the potential are given by 



V{z) = h 



a hgz + y^ j a k log(z 



Zk 



k=l 



(3.6) 



with h = —2ig. We treat h as a scaling parameter which relates the parameters of 
the Liouville theory with those of the matrix model. When /3 — 1, the integral (3.5) 
reduces to the Penner-type matrix model, where A/ is regarded as the eigenvalue of a 
hermitian matrix. The matrix model integral is not well defined unless the integration 
range and the Liouville parameter b are to be appropriately assigned. To fix this 
problem, one may consider the integrals with analytical continuation and put b 2 = —f3 2 
(or b = z\//3) so that integration is well defined even when any two of the integration 
variables coincide. This way one has the Penner matrix model (3.5) where iV is the 
size of the matrix determined by the neutrality condition (3.2). 

We now see how the stress tensor insertion in the conformal block 

(VooVo---V„) 

can be calculated in the Penner type matrix model, where (VqoVo 
conformal block associated with the correlator (3.1). For that, we first recall that the 
loop equation of the matrix model is obtained by changing variables as A/ — > Xi + 
and collecting terms of 0(e) in the partition function (3.5): 

N i n \ N 



(3.7) 

V n ) denotes the 



Aj — z 







E 

i,j=i 



(3 



(A; 



z)(Xj-z) 
1 



E 



7=1 



N 



1=1 



z — Xi 



+ 



? E 

i=i 



1-/3 
(A/ - zf 

V'(z) - V'iXj) 



z — Xi 



(3.8) 
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The first and second terms come from the variations of the measure and the Vander- 
monde determinant while the third and fourth terms come from the variation of the 
potential. When we define 

W(z ir -- ,z a ) = P ' 
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N 



y.- ! —-y.— 



V'(z) - V'jh) 

z — \j 



(3.9) 
(3.10) 



the equation (3.8) can be written as 



= g 2 W(z, z) + W{zf + V'(z)W(z) +g(V(3 



W'(z) 



(3.11) 



which is called the loop equation. For the Penner type model, the explicit form of the 
potential (3.6) implies that 



/(*) = V£ 



k=0 



z - z k 



d 

dk = log Z/j-Penner- 

dz k 



(3.12) 



Now, let us rewrite the loop equation (3.11) as 

V' 2 + f + hQV" = x{z) 2 + hQx'(z) - h 2 W{z, z) 



(3.13) 



where x(z) = 2W(z) + V'(z) and Q = b + 1/6. The left hand side of the equation is 
illuminating if one writes it as 



<p{z) 



where Z, 



cff 



h 2 

Z. 



/3— Pcnncr 

n < a <b<n( z a ~ z b) aaab • Then, equation ( |2.4[ ) and Z eS 
(T(z)VooVo---V n ) 



(VooVo • • • V n ) imply 

tp(z) 

In general, <p(z) has the form 



<p(z) 



(VooV ---V n ) 



P2n{z) 
Ylk=o( Z ~ Z k) 



(3.15) 



(3.16) 



where P^niz) is a polynomial of z with order In. This consideration allows us to 
consider the loop equation as the one generating the Virasoro constraint. The identi- 
fication (3.15) is consistent with the AGT relation Jl] as will be seen in section [5} 
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3.2 Colliding limit of the Penner models 



We now take the limit of a k — > oo, Zk — > in the matrix model (3.5) while keeping 

•an k 
r=0 r=0 0<si<---<s k <n (sj^r) i=l 

finite. In the CFT side, this gives a conformal block with an irregular vertex operator 
at z = and a regular vertex operator at z = oo, or equivalently, the inner product 
(Ri\I^). In the matrix model side, the same colliding limit gives the matrix model 
potential 

1 n 

- h V(z) = a\o g z-J2^ (3-18) 

k=l 

and the neutrality condition a + + bN = Q. We denote by Z M the partition 
function of the matrix model, that is, Z M = lim^^o.ar-^oo ^s-Penner- We call this type 
of matrix model "-D2n type" for a reason which we will see in section |5j This matrix 
model gives 

/W=-*'E^ (3-19) 

fc=0 

where v k = J22=i ac a+k-£r using the notation C£ = if £ > n+1. The term proportional 
to 1/z vanishes due to the identity V'(Xi)) = 0. 

Since we have taken the same colliding limit in the CFT side and in the matrix 
model side, we now identify 

Z M = (Ri\I (n) ), (3.20) 

where the right hand side is the inner product of an irregular vector \I^) of rank n 
and a regular vector \Ri). Although we have identified the regular conformal block 
with the rescaled partition function Z e s = Z^.penner Ylo<a<b<n( z a ~ Zb)~ 2oiaab , we here 
identify the irregular conformal block with Zm itself. This difference comes from the 
definition = lim^o.cv-j-oo {Yl < r<s < n (z r ~ z s) 2ar0la \Rn)) of the irregular vector. 

It is easy to note that for the potential (3.18) 

2n . 
k=0 

where A k = (k + l)Qck — Y^£=o c e c k-e with A = A a = a(Q — a){^] At the colliding 
limit, we obtain 

„ W . -1 (V- + / + hQV") = ± £ + g *S*M . (3.22) 

fc=0 fc=0 



^Recall here that cq = a. 
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Since <p(z) is identified with the stress tensor insertion into the conformal block, this 
shows that the non-vanishing expectation values of the Virasoro generators are read 
off as 

(Ri\L k \I^) J A fc when n < k < 2n 



(RtllW) I A fc + v k {\ogZ M ) when < k < n - 1 



(3.23) 



This is in perfect agreement with (2.6) and (2.8), and supports our identification (3.20). 



Thus, at the colliding limit the Penner model realizes the simultaneous eigenstate of 
the Virasoro generators L k with n < k < In. On the other hand, L k with k < n is 
represented in terms of differential operator v k . Thus, we will call the colliding limit 
of the Penner model as the irregular matrix model. 

3.3 Matrix models for general irregular conformal blocks 

The matrix model we have obtained above gives a Liouville conformal block with an 
irregular vertex operator at z = and/or a regular vertex operator at z = oo, or 
equivalently, an inner product (-Ri|J (n) ). The irregular vector is characterized 

by a, ci, • • • , c n and v fc(log Zm) for k — 0, • • • , n — 1, while the regular vector \Ri) is 
specified by aoo- Although the parameter does not appear in the matrix model 



potential (3.18), it is encoded in the matrix size A" through the condition a+a^+bN = 
Q. In particular, if = then we have no vertex operator insertion at z = oo. As 
will be seen in section [5j the matrix models for the potential (3.18) are related to the 



Argyres-Douglas theories of A 2n -3 and D 2n types. 

It is worth noting that the simultaneous eigenvalues of L 2 „, • • • , L n for the irregular 
vector |/( n )) are now encoded in the parameters in the matrix model potential V(z). 
On the other hand, the actions of L n _i, ■ ■ ■ , L on 1 1^) are specified by v k log(log Z M ) 



in the matrix model side. From (3.19), we find that ffc(log Zm) are determined when 
we fix f(z). Since fixing f(z) is equivalent to fixing the filling fractions of the matrix 
model, the distributions of A^ eigenvalues specify the actions of L n _i, • • • ,L on the 
irregular vector. 



By generalizing the argument in subsection 3^, we can now construct matrix mod- 
els for general Liouville conformal blocks with many regular and irregular vertex op- 
erators inserted. In fact, when we consider some additional regular /irregular vertex 
operators at z = z k for k = 1, - ■ ■ ,r, the matrix model potential is now given by 

i / n \ r / n k (k) 

h 




V(z) = [ a log z - £ -J-. ) + £ ( log(z -**)"£ ) > (3-24) 



with a neutrality condition a + Y7k=i a + a °° + bN = Q. The partition function 
of this matrix model reproduces the Liouville conformal block with (r + 2) vertex 
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operators inserted at z = 0, oo and z = z k . The vertex operator at z = z% is regular if 



(fc) 



for all j, while it is irregular if some c^' 0. 



It is also possible to make the vertex operator at z = oo irregular. One way to 
do this is to consider a colliding limit of some regular vertex operators into z = oo. 
However, in the matrix model side, we can easily see the effect of an irregular vertex 
operator at z = oo, just by changing the variables as Aj — > 1/Aj in the matrix model 
with (3.18). Under this transformation, the integrand of the matrix model integral 

N 



.1=1 



Ajf exp 



VP 

9 



(3.25) 



gives some changes. By exponentiating all the changes coming from the integration 
measures and the Vandermonde determinant, we obtain the same form of the matrix 
integral with a different potential 

1 



h 



V(z) 



log z - 



n u 
C k Z K 



k=l 



k 



(3.26) 



where the logarithmic term expresses the regular vertex operator now at z = and the 
other terms (together with oto) characterizes the irregular vertex operator at z — oo. 
Thus, in general, adding some polynomials of z to the potential leads to an irregular 
vertex operator at z — oo. For example, the matrix model for a Liouville conformal 
block with irregular vertex operators of rank m at z = and rank n at z = oo is given 

by 



V{z) 



a 



<°>los 




n J°°) 



k=l 



Z 



(3.27) 



with a neutrality condition + + bN = 
irregular singularity at origin, while a (oo) , 4 oo) 



- Q. Here a^\c^ characterizes the 
characterizes the one at infinity. Of 



course, we can further add some regular /irregular vertex operators at z = z k , just by 
adding the second bracket of (3.24) and modify the neutrality condition. 

Hence, if the matrix model potential V(z) is written as a sum of logarithmic and/or 
rational functions of z, the matrix model integral (3.25) gives a Liouville conformal 
block with regular and/or irregular vertex operators inserted. Note here that, among 
the large class of matrix models we have constructed here, a matrix model for two 
irregular states of rank 1 as well as a model for two regular states and one irregular 
state of rank 1 were already obtained in |llj . 



to d 



As will be seen in section |5| the matrix models for (3.24) and (3.27) are related 
4, M = 2 gauge theories which is not necessarily conformal. In particular, the 
potential (3.27) is associated with an asymptotically free theory denoted by v4 2m 2 n 
theory in the notation of [3]. 
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4 Irregular partition function 



The matrix models we have constructed above can be analyzed by the usual loop 
equation method with the genus expansion of /3-ensembles [28j[29j|30j. As an example, 



in subsections 4.1 and 4.2, we illustrate how to calculate the partition function of the 
D 2 n-type matrix models with potential (3.18), order by order in the genus expansion. 



In particular, we explicitly evaluate the first three expansion coefficients of log Zm for 
the -D4 matrix model. The result is in perfect agreement with the proposed ansatz 



(2.10) in the Liouville theory side. In subsection 4.3, we consider more general matrix 



models and argue that the small c 2 limit of the matrix models is consistent with (2.10) 
at the leading order of the C2-expansion. 



4.1 Property of the Dm type matrix model 

In this and next subsections, we concentrate on the D2 n matrix model with the poten- 



tial (3.18), which gives the inner product of a regular and an irregular conformal block 
(Ri\I^). We here denote by Z^ the partition function of the D 2n matrix model. 



One can evaluate the partition function using the loop equation 

^ = *W(z, z) + W{zf + V(z)W(z) + ^W'(z) 



(4.1) 



where f(z)/4 is now written as 
/(*) 



4 



h\ ^^(logzy, 



z 



2+k 



V k 



^ ac a+k 



dc a 



(4.2) 



k=0 o=l 

Recall here that we set c^. = for k > n. We can regard the loop equation (4.1) as a 



series of differential equations for Z M . To see this, we expand the quantities around 

z — 00 as 



W(z) 



hb 



g 2 W(z,z) 



£1=0 

hb 



E 

£,m>l 



E 

N 



Ef =1 (A/)0 



(4.3) 
(4.4) 



By collecting terms at each order of z in the loop equation (4.1), we obtain a system 
of differential equations. 

For example, in the case of n — 1 where the potential V(z) has two coupling 
constants a and ci, the loop equation gives us a single differential equation 



- v (log Z 



M ) 



A, 



(4.5) 
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where Vq = C\-M- and A 2 = (b 2 N 2 + 2baN — bQN). This differential equation solves 



dci 

the partition function as 



logi?j/ = — A 2 logci + constant. (4.6) 



where constant term is independent of a, ci,g,b and iV. 

In the case n — 2, the potential '/(z) includes three couplings a, Ci and c 2 , and the 
loop equation gives two differential equations 

-v Q.ogZ$)=A 2 (4.7) 

N 



-v 1 (logZ$)=2bNc 1 + A 3 (4. 



1=1 

where vq = c i^r + 2c 2 ^, t>i = c 2 ^, A 3 = 26 2 iV + 26a — 25Q, and A 2 is the same 
quantity as before. To find Z^ , one first notes that there is a homogeneous solution 
to vo = cit^j- + 2c 2 ^- since vo(H(t)) = for any function of £ = c 2 /c 2 . Thus, it is 

(2) A 

convenient to consider \ogZ M as a function of t,c 2 and solve (4.7) as 

\ogZ^ = -^\ogc 2 + H{t). (4.9) 
Applying V\ = c 2 ^ to log one has 

«x(logZg)) = -2 Cl t 2 ^|^ = _ 2ci t^'(t) . (4.10) 



Then, it follows from equation (4.8) that 



if'(t) = -L U 3 A/) + 2iV6 ) . (4.11) 



This shows that (X)j=i ^i)/ c i is the function of t only. Integrating over t one has the 
solution 

Nb A 3 f Gi(t) 



= -t + f y dt -ir (412) 

where = (X/jLi ^i)/ c i- Thus, the partition function is given in terms Gi(t). 

For n > 3 one finds n differential equations for Vk{log Zm) for < k < n — 1 
involving ((TrM a )) as well as ((TrM a ) (TrM 6 )) conn . Among them, we always have 

-v Q (\ogZ$)=A 2 . (4.13) 

Since any functions of the ratios 
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are in the kernel of Vq, the differential equation (4.13) implies 



log Z£> = log c„ + tf(|- ..,|). (4.15) 
The function H is determined by other differential equations. 



4.2 Details on D4 type 

In this subsection, we study the first non-trivial example of the n = 2 case in detail. 
The potential V(z) and the quantity f(z) are now given by 

\ V {z) = «togz-*-^ = * + * (4-16) 

n, 2; 2^ 2T" 5 



where we set d = — H 2 vo(logZ^) and di = — h 2 v\(\og Z$). As seen above, the 

7 (2) 

"* di dt' 



(2) 

partition function Z\J is given by 



iog4; ) (x,t) = -^io g c 2 + /d(t), ff(*) = J 



2h 2 d t' 2 



(4.17) 



Therefore, when we obtain d\j1c\ as a function of £ = c^jc 2 ^ we can evaluate the 

(2) 

partition function Z M . In fact, the quantity d\ is fixed by the filling fraction. Since 
the -D4 matrix model has two cuts, there is a single independent filling fraction. 

0(Fr 2 ) contribution 

To be more specific, let us consider the large N expansion of our matrix model while 
keeping h~N = 0(1), and evaluate \ogZ M order by order in h. In order to use the 
usual genus expansion method, we set a,Ci,c 2 = 0{N) so that V(z) is of 0(1). In 
other words, we rescale the parameters as 

« = -, c 1 = -, c 2 = -, (4.18) 

and treat a, ci,c~2 as of 0(1). In this setup, the resolvents and the partition function 
have the following expansions 



00 



W(z)= y %2W n (z)h n , W(z,z) = Y,W n (z,z)h n , \ogZ$ = Y J Fnh n ~ 2 . (4.19) 

n=0 n=0 n=0 

We can also expand do and d\ as 

00 

do = d 0;0 + fafoji, di = d ^ n , (4.20) 

n=0 
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with G?o ; o — bhN(bhN + 2a), g?o ; i — —bQhN. Then, the loop equation implies that at 
the leading order 



2W (z) + V'(z) = V^P-, (4.21) 

z 6 

Va(z) = (d , + a 2 )z A + (o?i ;0 + 2ac 1 )z 3 + {c\ + 2ac 2 )z 2 + 2c x c 2 z + c\. (4.22) 

The polynomial V^{z) has four roots, which implies that Wq(z) has two branch cuts 
on z-plane. Note here that two of the four roots of V&{z) is proportional to c 2 while 
the others are not. Let us denote a branch cut between the first two roots by T 2 , and 
a cut between the other two roots by IY Then, in the limit of c 2 —> 0, the branch cut 
T 2 shrinks into a point while I\ still has a finite width in z-plane. Then, we set the 
filling fraction condition as 

' 1 W(z)dz = (4.23) 



2m J Al 2 

where A\ is a cycle encircling the branch cut IY Since the total number N of eigen- 
values are fixed, we have a constraint N± + N 2 = N for 

^rjf W(z)dz=™N a , (4.24) 

where A 2 is a cycle encircling T 2 . Below, we see how the filling fraction condition 
(4.23) fixes the 0(h~ 2 ) contribution to the function H(t). 

We first assume Ni and N 2 are of O(N), which implies that the right hand side of 



(4.23) is of Oil). Then, the leading contribution to (4.23) gives 



hbNx = —£ W {z)dz = — I ^ V ^ dz. (4.25) 

7™ JAi 2711 J A! Z 

Since A\ encircles the cut I\ which does not shrink in the limit c 2 — > 0, we can expand 
the integral around c 2 = as 



Ai z3 JAx z3 

^ ^ 2ac 2 z 2 + 2c\c 2 z + (? 2 (2ac 2 z 2 + 2c\c 2 z + & 2 ) 2 



2QAz) 8QJz 



,2 



^^-dz + c 2 I -^±^dz + O(cj) (4.26) 



where 



Ai z3 J A! z 2 ^ Q A (z 



Q 4 (z) = (rf 0;0 + a 2 )z A + (di ;0 + 2ac x )z A + c{z l (4.27) 
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is independent of c 2 . Since a/ Qa(z) has only one branch cut r 1; we can evaluate the 



right-hand side of (4.26) just by taking residues at z = and z = oo, that is, 



1 t ^\dz = (bhN + ~ 3dl f 4 ~ 8& ' ld ^ c 2 + Q(% 



2m J Al z* V 2ci ; 8c? v 2; 

(4.28) 



From this and the filling fraction condition (4.25), we find 

= &£iV 2 + i|6aZVi(6^JVi + 2d) - 2bhN 2 {bhN 2 + a) + 26 2 fr 2 iV 1 iV 2 )t + £>(P), 

(4.29) 



where £ = c 2 /c 2 . Then, (4.17) implies that the 0(h 2 ) contribution to \ogZ^ is 

evaluated as 

fc _ 2r bN 2 bN(bN + 2a)^ 

h F = logc 2 

+ ^ {bNxibN! + 2a) - 2bN 2 (bN 2 + a) + 26 2 iV 1 iV 2 } log* + C(t),(4.30) 

up to constant. 
0{fr 1 ) contribution 

( 2 \ 

Let us evaluate the next-to-leading contribution to \ogZ y M '. By collecting terms of 
0{h) in the loop equation, we obtain 



1 (QzV(z) + zd 0>1 + d ljl ) - %\ogV A {z))' + (4.31) 



' ' 8^°'*™ 4z 

Then its period is evaluated as 

— / WAz)dz = — I — % (Qz 3 V"(z) + zdo-i + rfi-i) - -• (4.32) 

Here, the integral in the right-hand side can be expanded in powers of c 2 as 

dz 

' (Qz 3 V"(z)+zd 0]1 + d 1 . A ) (4.33) 



(do., - Qa)z + (d 1;1 - 2Qci) dz + &2 I ( 3Q 



dz 



Ai ^VQa(z) ' Ja 1 V ±z\[QAz) 

&(do-i — Qa)z 3 + (adi-i + cid ;i — 3<5«Ci)z 2 + ci(g? 1;1 — 2Qci),2 

4Q 4 W 3/2 

+0(c 2 2 ) (4.34) 
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The right-hand side of this equation is evaluated by picking up residues at z = and 



z = oo. Then, (|4.32|) is written as 
1 



2ni 



d 1:1 2c?d 0; i - 3di ;0 di ; i + 3Qcidi i0 - 4ciad 1; i 
Wl{z)dZ = ~W 1 + 86? " 



c 2 + 0(%) 

(4.35) 



Now, recall that we have the filling fraction condition (4.23). The 0(h ) contri- 



bution to (4.23) gives 



2m 



W x {z)dz = 0, 



(4.36) 



which implies that 



Qbh(2N 2 -N 1 ), - 



2ci 



-* + o(r). 



(4.37) 



Here we used d ;i — —bQhN and <ii :0 = 2c\bfiN 2 + 0(t). Then, by using (4.17), the 
(^(fir 1 ) contribution to logZj^ is evaluated as 



K ^i = ^— logc 2 + log* + £>(£), 



(4.38) 



up to constants. 



O(h ) contribution 

We now evaluate the O(ftP) correction to logZjj^ . The loop equation (4.1) at O(h ) 
tells us 



1,1:2 W (z,z) + 2W (z)W 2 (z) + W 1 (z) 2 + V'(z)W 2 (z) + ^-W[(z), (4.39) 



4z 3 

which implies that 

Wo(z) : 



2 



(4.40) 



Therefore, we first need to evaluate Wq(z,z) to calculate W 2 (z). 

In fact, Wq(z, z) is obtained from an another identity for the resolvents. By chang- 
ing the variables as 



N 



(4.41) 
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in the matrix integral (3.25) and collecting terms of 0(e), we obtain an identity 
= 



^W( Zl , Zl ,z 2 ) + (2W(z 1 ) + V'( Zl ))W( Zl ,z 2 ) + d W( Zl )-W( Z2 ) 

4 OZ 2 Z\ — Z 2 



hO 8 

-U( Zl ,z 2 ) + ^ — W(z u z 2 ), 



2 d Zl 



where 



V'( Zl ) - V'iXj) 



zi - A/ 

From the leading order of this identity, it follows that 

U Q (z,z)-Wg(z)/2 



, z 2 -\j 

J I connected 



W n (z,z) 



2W (z) + Viz) 



Furthermore, by using another identities 

\ I J 

U (z, z) can be rewritten as 

c 2 dW(z) W(z) 2W'(z 



_dW(z) 
dci 

W'(z), 

W(z) + zW'(z). 



U(z,z) 



z 3 dci 



Therefore Wn(z,z) can be written as 



W (z,z) 



c 2 dW (z) W (z) 2W' {z) Wg(z) 



^V 4 (z) \z 3 dc 1 z 2 
Then, equations ( |4.40[ ),( |4.31[ ),( |4.21| ) and ( |4.49[ ) lead to 
di-,2 z(l + g) 



WoZ 



(4.42) 
(4.43) 

(4.44) 

(4.45) 
(4.46) 
(4.47) 



(4.48) 



(4.49) 



4^/VJzj ^V^z) 

z (— 3c 2 . — c 2 z(3ci + az) + z 2 (c 2 + 6diaz + 3z (d 1;0 + 2z (<i 0;0 + a 2 )))) 
+ 6V 4 {z) 3 / 2 

r 

,z{ - 6(3Qc 2 + z(— + 2Qci - rf 0; i2 + Q&z)) 2 



96P 4 (^ 3/; 



+2QV ^c 2 + 2c 2 a + 3«(di; + 2c x a) + 6^ 2 (d Q;0 + a 2 ) J I, 



(4.50) 



up to total derivative terms which does not contribute to the period. Then the period 
£^ W 2 (z)dz has the following C2-expansion: 



2tu 



W 2 (z)dz 



3bh(4N 2 - (AN 2 - Ni)Q 2 



1 3bhN 2 + 2a 



t 



Ai 



2 2 

3(6 2 ft 2 (13iV 2 2 - 4N±N 2 - 2JV?) + 46^(4A^ 2 - jVi)d + 4d 2 ) ^ \ d 1;2 



-r 



+o{i 3 ). 



2ci 
(4.51) 



Then the filling fraction condition (4.23) implies that 



d 1;2 36fr(4iV 2 - (4iV 2 - NJQ* 



2ci 4 
This and (4.17) give the O(h ) correction to logZjff as 
_ 36(4iV 2 - (4iV 2 - NJQ 2 ^ 



P + 0{P 



(4.52) 



t + 0(t 2 



(4.53) 



up to constants. 



Partition function 



(2) 

From (4.30), (4. 38) and (4.53), we find that the partition function Z M of the D4 matrix 
model is written as 



7(2) 

J M 



bN(bN+2a-Q) 



l) e C2 WA (4. 



54) 



in the genus-one approximation of the genus expansion. 

Recall that we have identified this partition function with the irregular conformal 



block (Ri\I^). Since satisfies (2.10) and \Ri) is independent of c 2 , the partition 



function (4.54) should reproduce the small c 2 behavior of (2.10). Let us now check 



this. In the small c 2 limit, the partition function (4.54) is proportional to 

exp( 



-biV 2 (|6Ar 2 +2a-|Q) , bN 2 c\. 



C 2 



(4.55) 



When we identify —bN 2 = a — a', this factor correctly matches (2.10). The identifica- 
tion 



bNo 



i\/f3N 2 = a — a' 



(4.56) 



is reasonable because —bN 2 is originally the momentum of an intermediate state in 
the Liouville (regular) conformal block which should be identified with a — a'. 
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4.3 Small C2 limit in general 



In this subsection, we consider more general matrix model whose potential is of the 
form 



z Iz 



(4.57) 



Here G(z) is written as a sum of logarithmic and/or rational functions which are 
regular at z = 0. In the Liouville side, the partition function of this matrix model 
gives a conformal block (G|/( 2 )(ci, c 2 , a)) with an irregular vertex operator of rank 
2 at z = as well as regular and/or irregular vertex operators inserted away from 
z — 0. The vertex operators away from origin are characterized by G(z). We will 
not explicitly evaluate the partition function of this matrix model. Instead, we here 
consider the small c 2 limit of this matrix model and compare it with the proposed 
ansatz (2.10) of the rank 2 irregular vector |/^), including the shift of the momentum 
a — > a'. 



The matrix model with potential (4.57) generally have multiple cuts, among which 



a single special cut shrinks into a point in the limit of C2 — > 0. When the matrix 
model has p cuts, the spectral curve is the double cover of a Riemann sphere with p 
square-root cuts on it. Here, we denote by iVj the number of eigenvalues distributed 
on the i-th cut, and fix them with a condition Ym=i = ^ ■ This implies that, we 
choose the integration contour C, of the Ni eigenvalues so that it passes through the 
i-th cut and does not pass through the other cuts. Without loss of generality, we set 
the p-th cut to be the shrinking cut in the limit c 2 — > and therefore N p eigenvalues 
are distributed on the shrinking cut. Then the partition function of the matrix model 
is written as 




where C = (Ci, 



,C p _i), N=(N 1 

N, 



A(A) ^ Z Np (X) exp 
■ ■ ■ , N p _i) and 

\N\ N p 



9 



i=i 



(4.5J 



I lip \ |iV| lip 

JCp \j=l J 1=1 J=l 



N„ 



9 



j=i 



.(4.59) 



Below, we evaluate the leading contribution in the limit c<i — > 0. In the small c<i 
limit, the ^-integral becomes singular because N p eigenvalues are distributed on the 
shrinking cut, while the A-integral remains regular because no eigenvalue is distributed 
on it. In order to study the leading singularity in Zn p (\) we first rescale as £/ — > C2^i, 
which leads to 



V{c 2 z) 



1 

c 2 



1 

'2^2 



£i 

2 



alog(c 2 z) + G(c 2 z). 



(4.60) 
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Since G(z) is regular at z = 0, G(c 2 z) approaches to a finite constant in the limit of 
c 2 — > 0. Therefore, we can omit it in calculating the leading contribution to Z^ p (X). 
Note here that, since the width of the shrinking cut was of order 0(c 2 ) before rescaling, 
it is of 0(1) in £-plane after the rescalingj^] Then the leading contribution to Zn p (X) 
is now written as 



Zn p (X) 



^^N p ) 2 +iQ^pNp-2ia^N p 



\N\ 



\2f3N p 



exp F 



(4.61) 



i=i 



where the first factor comes from the logarithmic term in (4.60) as well as rescaling d^j 
and A(£) 2/3 . The second factor is the leading contribution from Ylj J (A_f — c 2 Cj) 2/3 - The 
remaining contributions are included in expF which is written as (up to a constant 
prefactor) 



p I=l 



expF = / (IIOA(0 2/3 exp 
with the reduced potential 



ic 2 



2z 2 



+ ac 2 logz. 



(4.62) 



(4.63) 



Note here that the term ac 2 \ogz gives a finite contribution to even i n 

the limit c 2 — > 0. 



What we need to do next is to evaluate the leading contribution to (4.62). For 



that, we use the loop equation for the reduced matrix model integral (4.62) 

/(*) 



(ic 2 /2) 2 W(z, z) + W{zf + (c 2 /2)QW'(z) + W{z)V'{z) - 
where we treat ic 2 /2 as a matrix model coupling constant to define 



0. 



W(z u 



P 



ic 2 /2 
VP 



2-s I N p 



1 1 



N n 



,1=1 



z-ti 



(4.64) 



(4.65) 



(4.66) 



Here and in the rest of this subsection, (O) stands for the expectation value of O with 
the reduced matrix integral (4.62). From the explicit expression for V and the identity 
(Efii V'fc)) = °> we obtain /(*) = do/z 2 + h/z 3 , with 



(4.67) 



1 All the other cuts now shrink into £ = oo. 
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Since the asymptotic behaviors of the resolvents are given by W(z) = (ic 2 /2)\/pN p / 'z+ 
(ic 2 /2)v^(^f p ii)jz 2 + O {z- 3 ) , W(z, z) = {z~ A ), the loop equation implies that 



d = -cl(pN 2 p -2i^aN p + </0iQN p Y (4.68) 
d x = 2 Cl c 2 i\fpN p - c\ (2(3N P -2ai^f(3 + 2^iQ\ (4.69) 



Note that these are exact expressions without any approximation. The equations 



(4.68) and (4.69) imply that all the quantities of the reduced matrix model (4.62) 



including expF itself, are completely determined by (5^i=i £/)■ F° r example, the free 
energy F is obtained by solving 



C-2 



dF d + c\d\ 
dc 2 



2c| 



dF 



d\ 



-.2 ' 



(4.70) 



Especially in the small c 2 limit, the eigenvalues are localized at a value that extrem- 
izes the potential V(£j), that is, = —l/c\. Note here that the Coulomb repulsion 
due to A (£) 2 is subleading and all the N p eigenvalues take the same value — \jc\ in the 
limit c 2 — > 0. The reason for this is that we here keep N p finite which is different from 
the usual 't Hooft expansion of the matrix model. By substituting £/ = — l/ci + C(c 2 ), 



we obtain (Xlr=i£.r) = —N p /c\ + 0(c 2 ), and therefore the equation (4.69) implies 

(2/3 iV 2 + 2{Q - a)i^N p )cl 



d x = 2 Cl c 2 i\f/3N p + 



Cl 



Then (4.70) is now written as 

OF i^N p c 2 PNp , iQVPN p 



c 2 



dc 2 
dF 



c 2 + 2 
2i^f3N pCl 



+ 



+ 0(c 2 ), 



(2f3N 2 p + 2 (Q - a) i^N^j i + O(o 2 



c 2 \ '/ Ci 

The solution to these differential equations in the limit c 2 — )■ is given by 



(4.71) 

(4.72) 
(4.73) 



WPNA 

c 2 



(2f3N 2 p +2(Q-a) i^N^ log Cl + 



f3N 2 tQ^N p 



+ 



log c 2 , 



(4.74) 



up to constant independent of C\ and c 2 . This leading behavior of F and equations 



(4.58), (4.61) imply that the original matrix model behaves in the small c 2 limit as 



-i^pN p (^i^pN p +2a-^Q) -i^pNp(-2i^pN p -2a+2Q) 



exp 



jT A(A) 2/3 exp 



c 2 



| A' | 



i=i 



(4.75) 
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where V(z) is given by 



•2g 



V(z) = (a + i^fflN^ log^ - - + G(z) 



(4.76) 



Note here that this potential is independent of c 2 and the coefficient of log z is shifted 
by iy/j3N p . 

The remaining matrix integral (4.75) is in fact identified with an irregular confor- 
mal block with one irregular singularity of rank one at z = and irregular/regular 
singularities away from z — 0, the latter of which is characterized by the function 
G(z). Hence, (4.75) implies that in the small c 2 limit 



(G\I (2 \ Cl ,c 2 ,a)) 



When we identify 



exp 



—i 



x(G\I {1 \ Cl ,a + t^N p )}. 



f3N p 



a 



a 



c 2 
(4.77) 



(4.78) 



we find that (4.77) is in perfect agreement with the leading term of (2.10), including 
the shift of the momentum a — > a'. In fact, the identification (4.78) is the expected 
one because i\ff3N p is interpreted as the internal momentum of the regular Liouville 
conformal block before taking the colliding limit. Since equation (4.77) holds for 
general G(z) as long as G(z) is regular at z — 0, it supports the ansatz (2.10) proposed 
in [3]. We here emphasize that the result of this subsection is valid for all orders of 
the genus expansion. 



5 Interpretation in gauge theories 

It was pointed out [U |3j that the irregular conformal blocks of the Liouville theory 
should reproduce the Nekrasov partition function of Ar gyres- Douglas theories, which 
is based on the observation that the colliding limits in the gauge theory side give 
Argyres-Douglas theories and some asymptotically free Af = 2 theories involving the 
Argyres-Douglas theories [32} [33] . Since we have already constructed matrix models for 
the irregular conformal blocks, we here explore the relation between our matrix models 
for irregular conformal blocks and the Argyres-Douglas theories in four dimensions. 

5.1 AGT relation 

In this subsection, we first review the relation between the regular conformal blocks 
and the Nekrasov partition functions of SU(2) super conformal linear quivers p.J. Let 
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cn 7 



a. 



7r-l 



I 



2 



I 



<5> 



I 



Figure 1: Left: The Liouville conformal block with (r + 2) regular vertex operators 
inserted. Right: The quiver diagram of the Af = 2, SU (2) gauge theory corresponding 
the left conformal block. Each circle stands for an SU(2) gauge group, and each box 
stands for an SU(2) flavor symmetry associated with a hypermultiplet. The external 
momenta of the Liouville vertex operators are now encoded in the mass parameters of 
the hypermultiplets, while the internal momenta are related to the Coulomb branch 
parameters associated with the gauge groups. 



us consider a Liouville conformal block 

Fao^'-^aJilkjAZk}) (5.1) 

with (r + 2) regular vertex operators V ak = : e 2ak ^ : inserted. Here {jk} specify the 
intermediate channels as in the left picture of figure [l] while {z k } denotes the loci of the 
vertex operator insertions. The AGT relation states that this conformal block gives a 
Nekrasov partition function 

Z Nek ( y m i ,a i ;e 1 ,e 2 ) (5.2) 

of a J\f = 2, SU (2) quiver gauge theory whose quiver diagram is given by figure [TJ 
Here the f2-background parameters e\ and €2 are related to the scale parameter H and 
the Liouville charge Q = b + 1/6 through 

e x = hb, e 2 = h/b. (5.3) 

The mass parameters TOj for the hypermultiplets are related to the external momenta 
in Liouville theory by 

m Q rrioo Q 

«o = -r- + tt> a oo = -r- + 77; oi k = — for k = l,---,r. (5.4) 

The internal momenta 7^ are related to the Coulomb branch parameters a k by 

7fc = ¥ + f f ° r fc = 1 ''"' r_1 - ( 5 - 5 ) 

The (r — 1) independent loci of the Vertex operator insertions are encoded in the (r — 1) 
gauge couplings e 2mTk of the quiver gauge theory. 
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The Seiberg-Witten curve of the SU(2) superconformal linear quiver theory is 
generally written as x 2 = 02 (z) where z is a coordinate on the Riemann sphere and 
the Seiberg-Witten differential is written as Asw = %dz. This 02 {z) is identified in p] 
with the stress tensor insertion into the conformal block 

(T(z)V ao (z )---V aoc ( Zoo )) 



lim — h 2 - 



where 



(V QO (^o) 

) stands for the conformal block. 



V aoo (Zoc)) 



(5.6) 



Here, we can see that the equation (3.15) in the matrix model is consistent with 



the identification (5.6). In fact, (3.15) means that the right-hand side of (5.6) is given 



by V'(z) 2 + f(z), and therefore the Seiberg-Witten curve is given by 

x 2 = V'(z) 2 + f(z), 



(5.7) 



in the matrix model side. The fact that ( |5.7 ) reproduces the correct Seiberg-Witten 
curve was already seen in [5J [EE] . 



Since the Penner type matrix models (3.5 ) reproduces the Liouville (regular) confor- 
mal blocks, it is expected that the partition function of the Penner type matrix model 
reproduces the corresponding Nekrasov partition function. The parameter identifica- 



tions are given by (5.3), (5.4), and the Coulomb branch parameters a& are identified 
as [51 [23] 

1 



a k 



W{z)dz 



HQ 

2 : 



(5.J 



where Ak is an appropriate k-th A-cycles of the spectral curve. 



5.2 Hitchin system with irregular singularities 

To see the effect of the colliding limits in gauge theory side, we now briefly review the 
six-dimensional origin of the d — 4,jV = 2 gauge theories. A class of d — 4,JV = 2 
supersymmetric gauge theories, including the above mentioned SU(2) superconformal 
quivers, is obtained by compactifying the six-dimensional (2, 0) supersymmetric A\ 
theory on a Riemann surface C, with a partial topological twist [211 Here the 
(2, 0) Ai theory is the low energy effective theory on the stack of two M5-branes, and 
the topological twist leaves eight supercharges in four dimensions. We can introduce 
codimension two half-BPS defects on the six-dimensional theory keeping the four- 
dimensional M = 2 supersymmetry. Such defects are point-like on C and give some 
punctures on it. In this paper, we only consider the case where C is a (punctured) 
Riemann sphere. The Coulomb branch B of the four-dimensional theory is parameter- 
ized by the vacuum expectation values of some chiral operators, which are encoded in 
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a quadratic differential 02 on C. Here the chirality in four dimensions implies that 02 
is holomorphic on Cjj Then, the Seiberg-Witten curve of the four-dimensional theory 
is written as x 2 dz 2 = 4>2{z) with the Seiberg-Witten differential xdz. 

The identification of B with the space of quadratic differentials is easily understood 
when we compactify the theory further on S 1 of radius R, which gives a d — 3, AT — 4 
supersymmetric theory The moduli space A4 of the three-dimensional theory is a 
torus fibration over B, where the fiber directions locally parameterize the electro- 
magnetic Wilson loops along the S 1 . Here M. can also be viewed as a moduli space of 
five-dimensional supersymmetric Yang-Mills theory compactified on which is the 
space of solutions to the Hitchin equations [32] 

F + R 2 [tp,<p]=Q, d A ip = d A Jp = Q, (5.9) 

on C with some boundary conditions at punctures, modulo gauge transformations. 
Here F is a curvature of the gauge connection A of a Si7(2)-bundle V on C, the "Higgs 
field" (f is an (End K)-valued (1, 0)-form on C, and d A = dz(d z +A z ), d A = dz(d z +A z ). 
It is known that the dimension of the Hitchin moduli space M. is twice the dimension 
of B. The projection of A4 to B is given by picking a unique Casiminar tr((p 2 ) of 
ip, which is identified with the quadratic differential of the four-dimensional Seiberg- 
Witten theory: 2 = tr((^ 2 ). 

The Coulomb branch B of the four- dimensional theory is now identified with the 
space of quadratic differential tr(y? 2 ) with fixed boundary conditions at punctures. 
Suppose that we have a puncture at z = 0. By trivializing the bundle V — > C near the 
puncture, the Higgs field generally has the following boundary condition near z = 0: 

9 ~ dZ V^ 1 + ^ + " ' + ~z + regular ) ( 5 - 10 ) 

up to gauge equivalence, where t k take values in a Cartan subalgebra of sZ(2,C). If 
n = then the singularity is called the "regular singularity," while if n > then 
it is called the "irregular singularity." In order to make tr(ip 2 ) singlevalued, we can 
consider nGNorneN+1/2. However, in this paper, we only consider the case of 
n G N. The boundary conditions tj. specify the singular behavior of the Higgs field, 
and interpreted as coupling constants and masses of the four-dimensional theory. In 



particular the mass parameter is encoded in to- The boundary condition (5.10) implies 



5 Although this 4>2 comes from the vev of some chiral operator in six dimensions, it is not a 
scalar but a differential on C, due to the topological twist to realize J\f = 2 supersymmetry in four 
dimensions. 

6 To be more precise, M. is the space of five-dimensional BPS configurations which is invariant 
under the three-dimensional Poincare transformations. 
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that the meromorphic quadratic differential is expanded near z = as 



Note here that the coefficient of l/2 fe+2 for n < k < 2n are completely fixed by the 
boundary conditions while those for k < n depend also on the regular terms in 



(5.10). This is the same situation as (2.6) in the Liouville side and as (3.22) in the 



matrix model side. 

When the defects are regular, they generate hypermultiplets in four dimensions. 
Each such defect is characterized by a single parameter to which is interpreted as the 
mass parameter of the hypermultiplet. In particular, if the A\ Hitchin system on 
C = P 1 has only regular singularities, then the four-dimensional theory has a weak- 
coupled description as a SU(2) superconformal linear quiver gauge theory, as in the 
right picture of figure [Tj The partition function of this gauge theory corresponds to a 
regular conformal block of the Liouville theory, via the AGT relation. 

Now, let us consider the colliding limit of the Liouville conformal block which we 
have reviewed in section [2j In the Liouville side, it gives an irregular conformal block 



satisfying the Ward identity of the form (2.6). Through the identification (5.6), we 



find that such an irregular conformal block is realized by irregular singularities in the 
Hitchin system on C . In particular, (n + 1) parameters a, C\, ■ ■ ■ , c n of the irregular 
vector of rank n correspond to the boundary conditions £q, • • • ,t n of the Higgs field 
in the Hitchin systemj^] Such irregular singularities in the Hitchin system are known 
to give Argyres-Douglas theories as well as some asymptotically free M = 2 gauge 
theories in four dimensions [32j [2j [33] . Then, it was pointed out jU [3] that the irregular 
conformal blocks of the Liouville theory should reproduce the partition function of the 
Argyres-Douglas type theories. 

Since we have already constructed matrix models which reproduce the irregular 
conformal blocks of the Liouville theory, we can then provide the matrix model real- 
ization of the partition functions of the Argyres-Douglas type theories. In the next 
subsection we first review the Hitchin system which realizes the Argyres-Douglas the- 



ories, and in subsection 5.4 we specify the matrix models for them. 



5.3 Argyres-Douglas theories 

The Argyres-Douglas type theories were originally discovered by taking the IR limit 
at a special point on the Coulomb branch of Af = 2 gauge theories [3J1 ESI GEE EZl EB] • 

7 Since the Cartan subalgebra of s£(2,C) is (complex) one-dimensional, each tk is determined by 
one complex parameter. 
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At the point on the Coulomb branch, some mutually non-local BPS particles become 
masslessj^] which suggests the IR theory is superconformal [22]. Since the theories do 
not have a Lagrangian description, it is not easy to perform a detailed study on them. 

Recently, it was pointed out that the Hitchin system with irregular singularities 
realize some of the Ar gyres- Douglas type theories [321 IH El E3]- Below, we briefly 
review the section 4 of [33] to describe how the A\ Hitchin system realizes some class 
of Argyres- Douglas theories which will turn out to be related to our matrix models. 

We first note that, in order for the Seiberg-Witten differential xdz to have scaling 
dimension one, there is a constraint [x] + [z] = 1. Although the SU(2) superconformal 
linear quivers realized by Hitchin system with regular singularities have [x] = 1 and 
[z] =0, the Argyres-Douglas theories have different scaling dimensions with [z] ^ 0. 
Then the superconformal symmetry of the theories implies that there are at most two 
singularities on the Riemann sphere C, that is, one at z = and the other at z = oo 

A 2n -3 type Argyres-Douglas 

Let us first consider the case with a single irregular singularity of degree n at z = oo 
and no regular singularities. In this case, the meromorphic quadratic differential <p 2 {z) 
is written as 

tr(<^ 2 ) = z 2n ~ 2 + u 2 z 2n ~ 4 + u 3 z 2n ~ 5 + ■■■ + u 2n . 2 , (5.12) 

where we rescaled z, <p and shifted z so that the coefficients of z 2n ~ 2 and z 2n ~ 3 are 1 
and 0, respectively]^] Here and in the rest of this paper, we omit dz 2 in the right-hand 
side. Then the curve x 2 = tr(<£> 2 ) gives the Seiberg-Witten curve of the (Ai, A 2n _ 3 ) 
type Argyres-Douglas theories in the notation of [40, 33J. In this paper, we just call 
them A 2n -3 type Argyres-Douglas theories. 

The parameters Uk stand for coupling constants and the vacuum expectation values 
of the corresponding relevant operators of the Argyres-Douglas theories. In particular, 
the superconformal point is given by = 0. The variables z and x has dimensions 
[z] = 1/n, [x] = 1 — \ jn so that [x] + [z] = 1 and 2[x] = (2n — 2)[z\. This fixes the 
scaling dimensions of Uk as [uk] — k/n. 

If n G N, the dimension one parameter u n is a mass deformation parameter of the 
theory, while u 2 , ■ • ■ , «n-i of dimension less than one are coupling constants associated 

8 By "mutually non-local" charges, we denote electro-magnetic charges T% and T2 whose Dirac- 
Schwinger-Zwanziger product (ri,^) is non- vanishing. See also |39j for a recent study on the Higgs 
branch of the Argyres-Douglas type theories. 

9 In this rescaling, we keep the Seiberg-Witten differential xdz invariant. This is always possible if 
n > 0. 
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with some relevant operators. The remaining parameters u n+ i, • • • , U2 n -2 are identified 
with the vacuum expectation values of the relevant operators, which follows from 
[uk] + [«2n-fc] = 2. On the other hand, if n G N+ 1/2 then there is no mass parameters. 
The parameters u 2 , ■ ■ ■ , ^n-1/2 are now couplings of the theory, and w n+ i/ 2 , ■ ■ ■ , u 2n -2 
are the vev of relevant operators. In this paper, we only consider nGN cases. 

£*2n type Argyres-Douglas 

Now, let us consider the Hitchin system with an irregular singularity of degree n at 
z = and a regular singularity at z = 00. The meromorphic quadratic differential is 
now given by 

^ ) = + ^2n+l + ■ ■ ■ + — + ^T> ( 5 ' 13 ) 

where we rescaled 2, <£> so that the coefficient of the first term is one. We have no 
freedom to shift z because we now have two punctures on C. The curve x 2 = tr(<£> 2 ) is 
equivalent to the Seiberg-Witten curve of the (A 1 ,D 2n ) type Argyres-Douglas theories 
in the notation of [lOj [33], which we just call Dm type Argyres-Douglas theories. 

The scaling dimensions of z and x are now [z] = —l/n,[x] = 1 + 1/n, which implies 
[uk] = k/n and [m] = 1. In addition to parameters U2, ■ ■ ■ ,U2 n -2 which have the 
same interpretation as for the A 2n -3 type theories, we now have U\, U2 n -\ and m. The 
condition [ui] + [u2 n -i] = 2 implies that U\ is a coupling constant associated with a 
relevant operator whose vev is given by u 2n -\- The parameter m is an additional mass 
deformation parameter associated to the regular defect at z = 00. 

What is important here is that the regular puncture at z = 00 is associated with 
a SU(2) flavor symmetry. Then, we can perform "gauging" the diagonal SU(2) flavor 
symmetry of Z^n and D2 m theories, by introducing an additional SU(2) vector multi- 
plet. Such an operation is interpreted as cutting a hole at each regular puncture and 
gluing them with a tube [31], which results in a Riemann sphere with two irregular 
singularities at z = and z = 00. Since we introduced an additional vector multi- 
plet, the resulting gauge theory is not conformal but asymptotically free. This type of 
theories is called v4 2m ,2n theory in [HI 0]. 

^2m,2n theories 

Let us briefly review the Seiberg-Witten curve of the A 2mj 2n theory, following [I]. 
Suppose we have an irregular singularities of degree m and n at z = and z = 00, 
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respectively. Then the meromorphic quadratic differential is now written as 



2n-3 i \2^2n-2 

i 

(5.14) 



where we rescaled z, tp so that the coefficients of l/z 2m+2 and z 2n ~ 2 are the same. The 
Seiberg-Witten curve of the theory is given by x 2 = tr(<y? 2 ). From our construction, 
Ui, 4 4 4 , u^m-i are parameters of the D 2m theory while Ui, • • • , u 2n -i are those of the 
D 2n theory. The additional parameters A and u are the dynamical scale and the 
Coulomb branch parameter for the additional vector multiplet. 

Theories involving AD theories as building blocks 

By generalizing the above argument, we can consider the gauge theories associated with 
the Ai Hitchin system with many regular and irregular singularities on C = P 1 . Then, 
the quadratic differential has the following singular behavior near each singularity, say 

at z — Zk' 

u (k) u {k) w£i 

t] V } = {z _ z ° k) 2m k+ 2 + {z _ z ly mk+ l +■■■+ {z T Zk] 2 + • ■ ■ ' ( 5 ' 15 ) 

where denotes the dynamical scale. The parameter u^\u[ k \--- ,Um k are fixed 
by the boundary conditions of the Higgs field tp at the puncture. The corresponding 
four-dimensional gauge theory is generically not conformal. The Seiberg-Witten curve 
of the theory is given by x 2 = tr((p 2 ). The complex structure moduli space of the punc- 
tured Riemann sphere is identified with the space of marginal gauge couplings,^ These 
gauge theories are called "wild quiver gauge theories" in jl], and can be constructed 
from SU{2) vector multiplets, bifundamental and fundamental hypermultiplets, and 
D 2 n type Argyres-Douglas theories. For example, if there are three irregular singulari- 
ties and two regular singularities on C = P 1 , then the resulting four-dimensional theory 
is described by a quiver diagram depicted in figure |2j in an appropriate weak-coupling 
regime of marginal couplings. 

The dimension of the Coulomb branch of the resulting four-dimensional gauge 
theory is given by jl] 



3 + t + 2% + 



' 2mk — 1 



(5.16) 



10 Note that generically not all the gauge couplings are exactly marginal. 
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Figure 2: Left: A Hitchin system on C = P 1 with two regular and three irregular 
singularities. Right: A quiver diagram for the gauge theories associated with the 
Hitchin system. We here take an appropriate weak coupling limit of the marginal 
gauge couplings. 



where £ and i are the numbers of regular and irregular punctures, respectively. When 
all nik are integers, this reduces to 

i 

-3 + (£ + i) + Y, m k- (5-17) 

k=l 

5.4 Matrix models for AD-theories 

We now specify matrix models which realize the above mentioned gauge theories, 
following the conjecture that the Liouville irregular conformal blocks should reproduce 
the partition functions of the Ar gyres- Douglas type theories. Since our matrix models 
are associated with irregular vectors of integer ranks, we here only consider the Hitchin 
system with irregular singularities of integer degrees. 

In particular, we explicitly write down the matrix model potentials for the A2 n -3, 
D 2n Argyres-Douglas theories and A 2m)2 „ theories for m, n G N. In this section, 
we only consider the case of m, n G N. We will briefly comment on the other cases 
m, n G N + l/2 in section[6j In the matrix model side, the spectral curve coincides with 
the Seiberg-Witten curve of the gauge theories. The coupling constants of Argyres- 
Douglas theories are encoded in the matrix model potentials, while the vev of the 
relevant operators are encoded in the eigenvalue distributions of the matrix models. 
The partition function of the matrix models are conjectured to give the Nekrasov 
partition functions of the corresponding gauge theories. 

Matrix models for A 2n -3 AD-theories 

The A2n-3 Argyres-Douglas theory is realized by the A\ Hitchin system with one 
irregular singularity at z = oo. Such a Hitchin system is related to the one point 
function of an irregular vertex operator of the Liouville theory. We here only consider 
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the case nGN. In the matrix model side, the corresponding potential is given by 

i n k 

s"m = -£ip (5 ' 18) 

fc=l 

By rescaling and shifting the eigenvalues A/, we can set c n = 1 and c n _i = 0. Such a 
rescaling just gives some constant multiplication to the partition function. With the 



potential (5.18), we have 

n-2 

(ArJ \ = 

z — X i 
with 

N /n-2 



/(^4 g ^(E rW :r (A,) )-X>^ (5-19) 



1=1 \ l=k 

„2 _ x//U\2 



Then, the spectral curve of the matrix model x = V'(z) + f(z) is exactly the same 
form as (5.12). In particular, the coefficients U2, • • • ,u n -i in ( 5.12[ ) is completely 



determined by the (n — 2) coupling constants in the potential ( 5.18[ ) p] The mass 
deformation parameter u n of the AD theory is determined by 

d n -2 = -2h 2 bN. (5.21) 

On the other hand, the other coefficients u n , u n +i, • • • , «2n-2 depend on c? 2) • • • , d n -\. 
In fact, these (n — 2) quantities are determined by (n — 2) filling fractions of the matrix 
model. The matrix model originally have iV eigenvalues, and they are distributed along 
the cuts of the spectral curve. Now, the matrix model has (n — 1) cuts, and therefore 
we have (n — 2) independent filling fractions: 



gy/PN k = I W(z)dz, (5.22) 



2ni 



where Ak is a cycle encircling the k-th cut and we have a constraint Ylk=i ^ k = 
Then, fixing all the filling fractions determines {dk}- 

Thus, we have found that the couplings and mass parameter of the AD theory is 
encoded in the couplings and the matrix size N of the matrix model. Since the matrix 
size N is related to the Liouville momenta through + bN = Q, this is consistent 
with the fact that the mass deformation parameter of the AD theory comes from to i n 



(5.10) in the Hitchin system. We have also found that the vacuum expectation values 
of the relevant operators in the AD theory are now encoded in the filling fractions of 
the matrix model. 



11 Recall that we are now setting c n = 1 and c„_i = 0. 
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We here briefly note that (5.18) is consistent with the original argument in 
In fact, the A 2n -3 type Ar gyres- Douglas theories can be geometrically engineered by 
type IIB string theory on a Calabi-Yau singularity 



uv + x 2 + z 2n ~ 2 = 0, 



(5.23) 



which has the matrix model realization with the potential V(z) = z n /n. Including the 
relevant deformations of the gauge theory now corresponds to deforming the potential 



as (5.18). In this paper, we instead derive (5.18) through the scaling limit of the 



Penner type matrix models which has a direct connection to the Nekrasov partition 
function of SU(2) superconformal linear quivers. 



Matrix models for D 2n AD-theories 

We now turn to the D2 n type Argyres-Douglas theories. The corresponding Hitchin 
system has an irregular singularity of degree n at z = and a regular singularity at 
z = oo. In the Liouville side, this setup corresponds to considering an inner product 
(.Ri|/ (n) ). We again concentrate on the case of n G N. Then, the corresponding matrix 
model has the potential (3.18), where we set c n = 1 by rescaling the eigenvalues A/. 



Then the quantity f(z) is now written as 



n-l 



m = £ 



yk+2 ' 



(5.24) 



fc=0 



with dk = — h 2 Vk{log Zm) in the notation of (3.19). The term proportional to 1/z 
vanishes due to the symmetry (J^jV'^Xj)) = 0. 

The spectral curve of the matrix model x 2 = V'(z) 2 + f(z) is now of the same form 
as (5.13). The couplings «i, • • • ,u n -i and mass parameter u n of the AD theory are 
completely fixed by n coupling constants in the matrix model potential (3.18l)F^ The 



mass parameter m 2 associated with the regular singularity is now determined by do. 
In fact, this do is fixed by the matrix size N. To see this, we recall the loop equation 

g 2 W{z, z) + W{z) 2 + ^-W'(z) + V'(7.)W(z) - K Z > 

Since W{z) = hbN/(2z) + 0(z~ 2 ), W(z, z) 
infinity, we obtain 



-W\z) + V'(z)W(z) 

0(z-% V'{z) 



4 



0. 



(5.25) 



ha/z + 0(z 2 ) around 



d = h 2 {b 2 N 2 - QbN + 2abN) 



(5.26) 



without any approximations. Thus, the mass parameter m in (5.13) is encoded in 

Q, this 



N in the matrix model. Since N is related to through a + + bN 
12 Recall that we are now setting c n = 1. 
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is consistent with the fact that m in (5.13) comes from the simple pole of the Higgs 
field in flOo| ). 

On the other hand, the vacuum expectation values of the relevant operators u n +i, ■ ■ ■ , «2n-i 
are determined by the remaining d k for k — 1, • • - , n — 1. The (n — 1) quantities are 
determined by (n — 1) filling fractions of the matrix model. In fact, since the spectral 
curve is of the form 



V 2n (z) 



- 2 = i^ (5-27) 

with some 2n-th order polynomial 7 ? 2n(^)j the matrix model has now n cuts. Therefore, 
there are (n — 1) independent filling fractions 



g^/3N k = ^-.i W(z)dz. (5.2* 



Thus, we have again found that the couplings and mass parameter of the AD theory 
is encoded in the couplings and the matrix size N of the matrix model, while the vev 
of the relevant operators in the AD theory are encoded in the filling fractions of the 
matrix model. 

Matrix models for A 2m ,2n theories 

The A 2m ^n theory is associated to the Hitchin system with two irregular singularities 
of degree m and n at z = and z = oo, respectively. We here assume m,n G N 
and m, n > 2. In the Liouville side, the partition function of this gauge theory is 
expected to give an inner product (/( m )|/( n ))|^j corresponding matrix model has 



the potential (13.27]), where we set c™ = = A by rescaling eigenvalues A/. Then 



the quantity f(z) is now written as 



m— 1 , 

/M = E % < 5 - 29 > 



where 



d k = -^ 2 bJ2\J2~hk) for -l<k<m-l, (5.30) 

N I n \ 

d k = -2h 2 bY" ( V c^Af * ) for - n < k < -2. (5.31) 



13 Note here that the ^2,2 theory is the SU(2) gauge theory with two flavors, and the matrix model 
realization for this theory was studied in . 
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The spectral curve of the matrix model x 2 = V'(z) 2 + f(z) is then exactly of 
the same form as (5.14). The couplings and mass parameters u\, • • • ,u m of the D 2m 
theory are completely fixed by the couplings a^°\ c[ \ ■ ■ ■ , c^Li of the matrix model, 
while the couplings U\, ■ ■ ■ , w„_i of the D 2n theory are fixed by c[ , ■ ■ ■ , in the 



potential (3.27). The mass deformation parameter u n of the D 2n theory is encoded in 
d_ n = —2h 2 AbN, that is, the matrix size N. The vev of the relevant operators in the 
D 2m and D 2n theories are respectively encoded in di, • • • ,d m -i and d-±, ■ ■ ■ , cf_ n+1 , 
while the Coulomb moduli u for the additional vector multiplet is determined by do. 
These (n + m — 1) moduli parameters are fixed by (n + m — 1) filling fractions of the 
matrix model. In fact, since the spectral curve is written in the form 

x- = (5-32) 

with some (2m + 2n)-th order polynomial V 2m +2n{z), the matrix model has now (m + 
n — 1) independent filling fractions: 



g^//3N k = ^- (b W{z)dz. (5.33) 

'Ah 



2m 



Matrix models for wild quivers 

When we consider an v^-Hitchm system with many regular and irregular singularities, 
the four- dimensional gauge theory is generically an asymptotically free theory involving 
Argyres-Douglas theories as building blocks. The Seiberg-Witten curve of the gauge 



theory is written in the form (5.15). The corresponding matrix model which gives a 



general irregular conformal block is described by the potential 

1 r ( m k S k ) \ m °o (OO) h 

where we assume Cm\ ^ if > 0. The singularity at z = z k is regular if = 
while it is irregular if > 0. Note here that we always have a (regular or irregular) 
singularity at z = oo unless the size of the matrix iV vanishes,^ 

The spectral curve x 2 = V'(z) 2 + f(z) of this matrix model is generally of the form 

^ = TT^ 7Z ^ \2m k +2> {5.35) 



where d = 2r + YTk=i ^ m k + 2m 00 — 2 and Vd{z) is a d-th order polynomial of z. This 
is easily shown when > 0. If = 0, we can see this as follows. First we note 



14 This is due to the fact that the resolvent is generally expanded around z = oo as W(z) 
hbN/(2z) + 0(z- 2 ). 
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that 



h 2 



Viz 



E 



Cm k 



+ 



H h 



d (k) 

u m k -l 



d 



(k) 



z - Z fc 

(*) 



-2 _|_..._|_ -1 



- Zfc) 



for some coefficients eft . Then, one might think that the spectral curve x 2 = V'(z) 2 + 
f(z) is of the form 



x 



V d+ i{z) 



(5.36) 



(5.37) 



nLi(^-^) 2mfc+2 ' 

where Vd+i{z) is a (2r + X^fc=i 2"T.fe — l)-th order polynomial of z. However, the 
coefficient of z d+1 in the numerator of (5.37) is 



(5.38) 



fc=i 



which turns out to vanish. In fact, this quantity is the residue of f(z) at infinity and 
equivalent to 



N 



-2hb(J2 v '( x i) ) =0 



(5.39) 



when rrioo = 0. Thus, the spectral curve is of the form (5.35) rather than (5.37) even 
when moo = 0. This fact implies that the number of independent filling fractions of 
the matrix model is always 



d 



- 1 



-3 + (r + l) + ^m fc + 



(5.40) 



k=l 



Since (r + 1) is the total number of (regular and irregular) singularities, 15 this is exactly 
the same as the dimension of the Coulomb branch of the corresponding gauge theory 



Partition functions of matrix models and gauge theories 

We have seen that our matrix models for irregular conformal blocks correctly reproduce 
the Seiberg-Witten curves of some Argyres-Douglas theories and wild quiver gauge 
theories. Recalling that the partition functions of the Penner type matrix models are 

15 Note that we always have a (regular or irregular) singularity at infinity. 
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conjectured to reproduce the Nekrasov partition functions of the SU (2) linear quivers, 
we now conjecture that the partition functions of our matrix models reproduce the 
Nekrasov partition functions of the corresponding Azn-z, -D2n-type Argyres-Douglas 



theories and wild quiver gauge theories. The parameters are identified as (5.3) and 



(0) mo Q M = moo Q (k) ™k , , 

h + 2 ' h r h' { 1 

where mo, m^, m k stand for some mass parameters of the gauge theory. The Coulomb 
branch parameters are identified as 

a k = -l W(z)dz + ^, (a D ) k = -l W(z)dz + ^, (5.42) 

where we take the 1-cycles of the spectral curve so that their intersections are given 
by 

(Aj, B k ) = S jk , (Aj, A k ) = (Bj, B k ) = 0. (5.43) 
In the genus expansion of our matrix model 

oo 

\ogZ M = Y,F n {ihy-\ (5.44) 

n=0 

the leading term F should particularly give the prepotential of the corresponding 
gauge theory. In fact, from the general property of the matrix model, it follows that 
Fq satisfies the special geometry relation: 

dF 



da k 



(a D ) k , (5.45) 



which is necessary in the gauge theory side. Thanks to this relations, Fq is determined 
by the spectral curve x 2 = 02 (z) and the meromorphic one-form xdz = (2Wo + V')dz. 
Since we have already checked that the spectral curve of our matrix models coincides 
with the Seiberg-Witten curve of the corresponding gauge theories, at least we can see 
that F correctly describes the IR physics of the corresponding gauge theories. It is 



worth studying the higher order terms of (5.44) further 



6 Summary and discussions 

In this paper, we have constructed matrix models which reproduce the irregular con- 
formal blocks of the Liouville theory on sphere. We have studied the matrix side of 
the colliding limit of the Liouville vertex operators, and pointed out that if the ma- 
trix model potential is written as a sum of logarithmic and/or rational functions then 



37 



its partition function reproduces a conformal block with insertions of regular and/or 
irregular states of the Liouville theory on sphere. In section [4| we have particularly 
studied the Z?2n-type matrix model in detail, and show that the partition function of 
the matrix model correctly reproduces the inner product of a regular and an irregular 
states. We have also shown that our matrix models generally reproduce the small c 2 
behavior of the irregular state \I^ 2 \ci, c 2 ; a)) proposed in [3J. In section [HJ we have also 
discussed the relation between our matrix models and the Argyres-Douglas theories 
in four dimensions. We have shown that our matrix models correctly reproduce the 
Seiberg-Witten curves of the corresponding gauge theories. 

We should here mention that we have not studied irregular singularities of half- 
integer degree in this paper. Especially, we have not studied matrix models for A 2n 
or D 2 n+i-type Argyres-Douglas theory. In fact, these theories cannot be realized by 
logarithmic or rational potentials of the matrix model. For example, the Seiberg- 
Witten curve of the A 2n Argyres-Douglas theory is of the form 

x 2 = z 2n+1 + u 2 z 2n ~ l + ■■■ + u 2n+1 , (6.1) 

and if this is realized as the (planar) spectral curve of the matrix model x 2 = V'(z) 2 + 
f(z) then it seems likely that the potential V(z) should have a square-root term z n+1 ^ 2 . 
However, it is not straightforward to generalize our argument to such a square-root 
potential. This complication comes from a different singular behavior near the irregular 
singularity of half-integer degree, which needs to be studied further. Note also that 
the colliding limits of Liouville vertex operators have not yet been well-established for 
irregular singularities of half-integer degrees. 

For future works, it would be interesting to study the higher orders of the c 2 - 
expansion (2.10 ) in the matrix model side, which will lead to a matrix model expression 
of an inner product (G\I^) with a generalized descendant \I%k), and to extend the 
method to the half-integer rank case. 

The application of our matrix models to the quantization problem of Hitchin system 
is also an interesting future problem. In [13] , the Penner type matrix models were used 
to quantize the Hitchin system with regular singularities. Since our matrix models can 
take into account irregular singularities in the Hitchin system, it would be interesting 
to generalize the argument in jl3] by using our matrix models. 

It is also worth studying the generalization to the higher rank of gauge groups. 
As pointed out in [5], the higher rank gauge groups correspond to the /3-ensemble of 
multi matrix models. By generalizing our argument to the multi matrix models, we 
can construct matrix models for A n Hitchin system with irregular singularities. Such a 
generalization will give a matrix model realization of (A n , v4 m )-type Argyres-Douglas 
theories in the notation of 140 1. 
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